Shortcuts to adiabaticity in a time-dependent box 
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A method is proposed to drive an ultrafast non-adiabatic dynamics of an ultracold gas trapped in a time- 
dependent box potential. The resulting state is free from spurious excitations associated with the breakdown 
of adiabaticity, and preserves the quantum correlations of the initial state up to a scaling factor. The process 
relies on the existence of an adiabatic invariant and the inversion of the dynamical self-similar scaling law 
dictated by it. Its physical implementation generally requires the use of an auxiliary expulsive potential. The 
method is extended to a broad family of interacting many-body systems. As illustrative examples we consider 
the ultrafast expansion of a Tonks-Girardeau gas and of Bose-Einstein condensates in different dimensions, 
where the method exhibits an excellent robustness against different regimes of interactions and the features of 
an experimentally realizable box potential. 



"Fast good" is a new culinary concept envisioned by the 
chef F. Adria aimed at creating a diet enjoying of two ap- 
parently mutually-exclusive features: fast-service and high- 
quality. When similar ideas are invoked in the quantum realm, 
one faces the adiabatic theorem which imposes a price. The 
preparation of a given target state with high-fidelity, free from 
spurious excitations, generally demands a long time of evo- 
lution and the implementation of an adiabatic dynamics. As 
a result, it comes as no surprise that the recent development 
of shortcuts to adiabaticity (STA) [ ], has led to a surge of 
theoretical [2-10] and experimental activity [11—13]. 

It is the purpose of this work to show that STA can be im- 
plemented in a time-dependent box, the paradigmatic model 
of a quantum piston [ ]. The relevance of this confinement 
is enhanced by the development of experimental techniques 
to create optical billiards for ultracold gases [15, 16], paint 
arbitrary potential shapes for Bose-Einstein condensates [17], 
and realize all-optical boxes [18] and analogous traps in atom 
chips [ ]. It has the potential to greatly advance the field of 
quantum simulation, facilitating the connection between ul- 
tracold atom experiments and condensed matter systems. 

The prospects of finding STA for box confinement seem 
challenging at the very least, based on the following consid- 
erations. Since the early insight by M. Moshinsky in 1952, 
it is known that the rapid expansion of matter-waves initially 
localized in a region of space exhibits quantum transients and 
excitations in the form of density ripples, a phenomenon re- 
ferred to as diffraction in time [20-22]. Moreover, reflections 
from the walls of the trap generally lead to Talbot oscillations 
and the formation of quantum carpet in the time evolution of 
the density profile [23, 24]. On top of that, a series of work 
during the last decades have shown that the suppression of 
excitations in an expanding box is inevitably constrained by 
the adiabatic theorem both in the non-interacting [25-28] and 
mean-field regime [29]. 

In spite of these results preventing ultrafast excitation-free 
dynamics in box traps, we show that fast non-adiabatic expan- 
sion or compression is in fact possible, allowing preparation 
in a predetermined finite time of the same target state than the 
adiabatic evolution. The key to this shortcut to adiabaticity is 
the use of an auxiliary potential term superimposed on the box 
trap. 



We shall start by introducing a dynamical invariant in a 
time-dependent box trap at the single-particle level, and use 
it to derive a scaling law that governs the time evolution. The 
inversion of this scaling law will be the key to engineer the 
STA. This result will be extended to many-body systems with 
a broad family of interactions. Finally, we shall illustrate the 
method with some examples relevant to experiments with ul- 
tracold gases, showing that the STA is robust in the presence 
of interactions and experimental imperfections. 
Dynamical invariants and self -similar dynamics 

For a time-dependent box of width ^(f), the scaling laws 
governing the dynamics of the expanding eigenstates reported 
to date are associated with trajectories of the form <i;(f) = 

[at 2 + bt + c] 2 (with a, b, c real constants) [25-28] , which 
turn out to be unsuited for engineering a STA (see below). 
Nonetheless, given a time-dependent Hamiltonian //(f), it is 
possible to build a dynamical invariant [30], J?(t ) such that 
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with spectral decomposition J^(t) = L«An|0 n (f))(0n(OI m 
terms of the set of eigenmodes |0„(f)) with eigenvalues X„. 
This is a particularly useful basis to describe the time evo- 
lution of an initial state *¥, by the superposition ^(x,?) = 
£„exp(/a„)|0„(f)), where the Lewis-Riesenfeld phase is 
given by a n = y Q {fy n (t')\ihd t i - H(t')\(j) n (t'))dt' fh, and can 
be understood as the sum of the dynamical phase and the 
Aharanov-Anandan phase [28]. For a time-dependent box of 
width ^ (t ) and initial width £, (0) = a dynamical invariant 
exists [31], 
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sin[nnx/%(t)] and eigenvalues 

A„ = ^yp, with k„ = n% /|o- The Lewis-Riesenfeld 
phase can be computed to be a„(t) = — n " \ 77(f) with 

77(f) = ^dt't^/t, 2 ^'). The condition in Eq. (1) for J to 
be an invariant requires the box potential to be supplemented 
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FIG. 1. Time-modulation of the potential trap along a shortcut to adi- 
abaticity. (A) Evolution of the box width t, (t) of a time-dependent 
box during a shortcut to different adiabatic expansions and a com- 
pression in a finite time T. (B) Frequency of the auxiliary exter- 
nal harmonic potential required to assist the self-similar dynamics, 
changing character from expulsive to attractive along the expansion 
process. The opposite sequence is required for a compression. The 
transient expulsive potential is responsible for the speed-up of the 
process. 



with an auxiliary harmonic term 
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whose frequency Q(t) = — is dictated by the ratio of 

the acceleration of the trajectory ^(t), and the trajectory it- 
self. Note that E,(t) > 0, so that for |(f) > 0, Q(f) is purely 
imaginary, and the auxiliary term t/ aux is a repulsive harmonic 
potential. If |(0) = 0, [H(0),J(0)] = so 6(0) and J(0) 
have common eigenstates. Further, if £, (0) = holds as well, 
then J/ aux (x,0) = and an eigenstate 0) of the box at 
t = evolves into v ¥„(x,t) = exp[/a„(?)]0 n (f), a key obser- 
vation to engineer a STA as we shall see. We note that the 
experimental implementation of <V aux (x,f) can be assisted by 
the same techniques used to create the box potential: the use 
of a blue-detuned laser [32] or direct painting of the required 
trap [15-17]. 

Shortcuts to adiabaticity: inverting the dynamical scaling 
law 

We next discuss how to implement a non-adiabatic expan- 
sion of the box by a factor /(t) = £, (t) /E,q in a given finite- 
time T suppressing excitations in the final state. We shall im- 
pose the condition t/ aux (x,0) = C/ aux (x, t) = 0. As in the adia- 
batic case, in a STA the time evolution of an eigenmode of the 
initial box should reproduce an eigenmode of the final trap. 
As at t = 0, this can be enforced by imposing the condition 
| (t) = (t) = 0. The set of boundary conditions at t = 0, T 
excludes the possibility of a linear ramp, as well as the fam- 
ily of trajectories, E,(t) — [at 2 + bt + c] 2 , considered so far 
in the literature [25-28]. However, it suffices to determine a 



polynomial ansatz for the trajectory ^(f) = Y.i ^, i. e. a 
scaling factor of the form y(f ) = = 1 + [f(x) - l]i 3 [10 + 
3s(2s — 5)], with s = t/i. This further determines the required 
time-dependent frequency of the auxiliary harmonic potential 
t/ aux (x,f) according to Eq. (3), 
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The trajectory, displayed in Fig. 1 shows that during an ex- 
pansion t/ aux (x,f) becomes an expulsive potential in an early 
stage (f < t/2), providing the speed-up required to achieve 
the STA in an arbitrary finite time T (bounds in the presence 
of perturbations will be discussed below). In a subsequent 
stage, t > t/2, £2 2 (f) changes sign and i/ aux (x,f) becomes a 
trapping potential, slowing down the expanding mode and re- 
ducing it to an eigenstate of the final Hamiltonian at t = T. 
Precisely the opposite behavior is exhibited during a fast non- 
adiabatic compression. Provided that an arbitrary £i 2 (f ) de- 
pendence can be implemented, a STA has no lower bound for 
T (notice however that Q.(t) ~ T _1 ). By contrast, the adiabatic 
condition 
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leads to the requirement mnkS, (t)t; (t) 

Note that the energy of the expanding mode 
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has two contributions, the first one being the adiabatic energy 
E n (t) = E n (0)% (0) 2 /£ (f ) 2 , and the second one depending ex- 
plicitly on ^, so that (H(t))„ = E n (x) given that a STA de- 
mands | (t) =0. This relation illustrates the fact that a STA is 
associated with a non-adiabatic evolution, which reproduces 
the adiabatic result at the end of the process. Moreover, STA 
work as well for excited states: the time evolution of the n- 
th eigenstate of the initial trap leads to the n-th eigenstate of 
the final trap at t = T. As a result, STA in boxes pave the 
way for fast population-preserving cooling in the following 
sense. Given a system described by the canonical ensemble 
with a density matrix e^^^/Trfe - ^ 11 ], where /3 = l/kgT, kg 
the Boltzmann constant and T the temperature, the final tem- 
perature reads 
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These results strictly hold for a box with infinite walls at 
x g [0,£(f)]. However, in the all-optical trap reported in [18], 
the end-cap lasers providing the box walls have a Gaussian 
profile which smooths the potential in a length scale a, i.e. a 
box trap of the form 
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A similar smoothing is expected in other physical realizations 
[17], Since this smoothing is expected to be the most signif- 
icant deviation of laboratory potentials from the ideal infinite 
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FIG. 2. Space-time density plot of the scaled density profile 
n(x, t)y(t) during a shortcut to an adiabatic expansion for the ground- 
state mode of a box with smooth boundaries. A high fidelity-state 
results from the STA (left) in a box created with Gaussian barriers 
of height V = I0 4 h 2 /m^ and a = | /20 (y = 10). In the absence 
of t/ aux (x) (center), and for a linear ramp box width t, (t) (right) the 
cloud is excited and exhibits fringes associated with the reflections 
of the boundary. The color coding varies from purple to red as the 
function value increases. The dashed line follows the trajectory of 
the moving wall ^ (t). 



box, it is important to consider its effect on the self-similar 
dynamics required for a STA. This can be quantified by the 
overlap between the states resulting from the expansion in an 
idealized and realistic box trap, |*P«(0) and |SP n< y(f)) respec- 
tively. Clearly the role of a decreases (increases) during a 
expansion (contraction) of the box. The numerical solution of 
the time-dependent Schrodinger equation for this box poten- 
tial is shown in Fig. 2, where the STA is compared with both 
the polynomial and linear expansion of the box in the absence 
of t/ aux (x,f). It is seen that the STA is the only successful 
strategy and that the process is robust even for a substantial 
smoothing of the potential barriers, where the target states de- 
viate from those of an idealized box. 

We have further explored numerically the dynamics under 
"concatenated STA", in which the overall expansion is split- 
ted into a sequence of k STA with either constant expansion 
factor 7 or constant box size increment between consecutive 
steps. The efficiency of the process exhibits a non-monotonic 
improvement with increaisng k, suggesting a natural scenario 
where to combine STA techniques with optimal quantum con- 
trol [8, 9]. 

Beyond adiabatic invariants: Shortcuts to adiabaticity in 
interacting many-body systems 

Knowledge of the adiabatic invariants for a single-particle 
in a time-dependent box has provided us with the insight of 
assisting the dynamics with an auxiliary potential to design 
STA for expansions and compressions in a finite time, with- 
out inducing quantum transients associated with diffraction in 
time, which are ubiquitous in this type of scenario [22]. The 
technique, can be directly applied to non-interacting gases and 
other many-body quantum fluids which can be mapped to non- 
interacting systems. It further suppresses the Talbot dynamics 
associated with quantum carpets woven by the density profile 
typically observed in boxes [23], and the question naturally 
arises as to its applicability to interacting systems [24]. The 
presence of interactions, e.g. a two-body potential, hinders 
the exploitation of the superposition principle in terms of the 
eigenmodes of the Lewis -Riesenfeld invariant. However, we 



note that to design a STA it suffices to enforce a self-similar 
dynamics and ultimately, no knowledge of adiabatic invariants 
is required. As a result, we next consider a broader family of 
many-body systems, confined in a box, defined by the Hamil- 
tonian 
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where q ; - <E MP, A qj is the Laplace operator in dimension D, 
the auxiliary term is now given by 
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and the two-body interaction potential obeys V(Aq) = 
A~"V(q), e.g. for the Fermi-Huang pseudopotential describ- 
ing s-wave scattering in ultracold gases, a = D. For a hard- 
wall box, r; = |q|; £ [0, we shall relax below this ap- 
proximation and consider realistic potential boxes as those 
created in all-optical setups . The case D = 1 corresponds 
to a box with one-wall moving (the symmetric case in which 
both walls move in opposite directions can be obtained by 
a Duru transformation [33]). For D = 2,3 cylindrical and 
spherical symmetry is assumed respectively. Without loss of 
generality, we choose the dimensionless time-dependent cou- 
pling constant £ = e(t) to satisfy e(0) = 1. A stationary state 
^(f) = ^(qi, . . . ,qw;0 of N particles and chemical potential 
ji follows for t > the evolution 
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with the boundary conditions *¥(t) = for |q|,- = %(t) (z = 
l,...,N, in addition to *P(f) =0for |q| ; = 0inD = 1), as long 

as 



e(0 = 7(0" 
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which can be implemented exploiting a Feshbach resonance 
or modulating the transverse confinement in anisotropic sys- 
tems [4, 34]. 

The self-similar dynamics in a STA leads to a scaling 
of all local correlation functions. In particular the den- 
sity of a given many-body state follows the law n(a,t) = 

Jc/q 2 ...rfqiv| , I , (q,q2,...,q7v;0| 2 = »[q/7(0,0]/y(0° By 
contrast, non-local correlation functions exhibit a non-trivial 
dynamics. The one-body reduced density matrix p\ (q,q';f ) = 
Njdq 2 ... dqn^iq, q 2 . . . , qjv; 0^* (q', q2 • • • , qw; of a state 
obeying Eq. (11), follows the scaling law pi(q,q';f) = 

exp [z m( ^2^ 2 ^ ] Pi (q/y, q'/7» °)» analogous to that ob- 
served under harmonic confinement [43]. The additional 
phase factor induces a major distortion of the momen- 
tum distribution, n(k,f) = Jc/qc/q'pi(q,q';f)exp[zk • (q — 
q')]/ (2n) D . However, a STA ensures that at t = % the Lewis- 
Riesenfeld phase factor vanishes, so that the final state ex- 
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FIG. 3. Time evolution of the momentum distribution of a Tonks- 
Girardeau gas confined in box during an ultrafast expansion and 
compression providing a shortcut to adiabaticity. In a shortcut to 
an adiabatic expansion (A), due to the auxiliary external potential 
in Eq. (3), the cloud is accelerated in the first half of the process. 
At s = 1 /2 there is a partial mapping of the momentum distribution 
to the density profile of the initial state of the cloud. In the second 
half of the process, the cloud is slowed down and at t = % it matches 
the result from an adiabatic dynamics: the state exhibits exactly the 
same correlations of the initial state, scaled up by the expansion fac- 
tor /(t) = 10 (N = 10). The reverse sequence is observed in a short- 
cut to an adiabatic compression (B), where the target state is reached 
with high-fidelity already for s = t jx = 0.75, after which the dynam- 
ics is actually adiabatic (/(t) = 1/10). 



hibits the same correlations of the initial state scaled by a fac- 
tor 7(t), 
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Examples 

In the following we shall illustrate different aspects of 
shortcuts to adiabaticity in some paradigmatic models. 

We shall first consider the evolution of correlations in a one- 
dimensional cloud of ultracold bosons in the limit of hard- 
core contact interactions, this is, in the Tonks-Girardeau (TG) 
regime [ 3 )]. This system, as well as its lattice-version, has 
become a favorite test-bed to study the breakdown of thermal- 
ization and adiabaticity [36]. Its many -body ground state is 



given by the Bose-Fermi mapping [35], 1 I , tg(?i , ■ ■ ■ ,Qn) 
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(-1) if x > (< 0) and e(0) =0. In a STA, the self- 
similar dynamics is inherited from the single-particle orbitals 
l ¥j(qi ( ,t) whence it follows that no tuning of interactions is 
required. Its density follows a scaling law for all t . The same 
holds true for the entanglement entropy with respect to a bi- 
partition [0, «<!;(?)] (a < 1) [37], The self-similar dynamics 
breaks down for the momentum distribution, which can be 
computed efficiently following [38], and we shall focus on 
its evolution along a STA. Different snapshots are depicted in 
Fig. 3A, and confirm that during an expansion the cloud is ac- 
celerated during the interval [0, t/2] and slowed down during 
[t/2,t]. The reverse sequence, is observed in a fast friction- 
less compression. The axis are scaled up by the expansion fac- 
tor y(f) in such a way that for an adiabatic dynamics, curves 
at different times would collapse into a single curve. Along 
a STA, the width and mean of the momentum distribution do 
not remain constant and change along the process. 

A similar distortion of correlations, known as dynamical 
fermionization, occurs in the dynamics of a cloud suddenly 
released from an arbitrary trap [39]. Under ballistic dynam- 
ics the asymptotic momentum distribution in a ID expansion 
evolves to that of the dual system, a spin polarized Fermi gas. 
In particular for a cloud released from a box the exact time- 
evolution is not self-similar [ ] but dynamical fermionization 
is observed [41, 42]. However, under a self-similar scaling 
law, the asymptotic n(k) maps to the density profile of the ini- 
tial state [43] and no dynamical fermionization occurs. This 
is the case of relevance to STA, where the dynamical scaling 
law in Eq. (11) holds. (We note that the case of the initial 
harmonic confinement is singular in that the free expansion 
is self-similar and that the the single-particle eigenstates can 
be written in terms of Hermite polynomials, which are eigen- 
functions of the continuous Fourier transform. As a result the 
asymptotic momentum distribution can be related to both the 
initial density profile and the momentum distrubtion of non- 
interacting fermions. See [5] for a discussion of STA in har- 
monic traps.) Moreover, this distortion of correlations is not 
restricted to expansion processes. Along a STA, this is shown 
in Fig. 3 for both expansions (A) and compressions (B). This 
is a spurious effect for the purpose of STA, which is to re- 
produce the adiabatic result in a finite short time. Indeed, the 
distortion induced during the first half of the STA associated 
with the accelerated expansion or compression, is compen- 
sated in the second half of the dynamics, in such a way that 
the correlations of the initial state are reconstructed at t = T 
and scaled by a factor 7(f). 

We next turn our attention to the design of STA for a BEC 
in time-dependent box trap, where different strategies can be 
adopted depending on the dimensionality and the regime of 
interactions. The time-dependent Gross-Pitaevski equation 
(TDGPE) governs the evolution of the normalized condensate 
wavefunction 4>(q, t ), 
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FIG. 4. Shortcut to an adiabatic expansion of a Bose-Einstein con- 
densate. The upper row shows the space-time evolution of a quasi- 
1D BEC (n(x,t)y(t)) while the lower row corresponds to a quasi-2D 
(«(r, f)y(f) 2 ) BEC cloud. The shortcut to adiabaticity (left) is com- 
pared with a polynomial (center) and a linear ramp (right) of the 
box width §(f) in the absence of the auxiliary harmonic potential 
(<t/£o = 1/20, g D = Wh 2 /m^ D , and y= 10). 



with |q| 6 [0,£(f)], for which adiabaticity conditions have 
been reported in [29]. The ansatz 



<J>(f) = 7 7 exp 



: m\q\ 2 y gqfr) 
2y7z ' ft 



4>[q/ 7 (0;0], 



(15) 



satisfies the TDGPE provided that 
dt' 



7(f) 



o 2 (A- ^ „ 

,, 2 - SD\f)--^D- 



(16) 



These relations constitute the box analogue of the well- 
known Castin-Dum-Kagan-Surkov-Shlyapnikov relations in 
harmonic traps [44, 45]. The two-dimensional case is special 
since the scaling law holds when g2D{t) is kept constant. 

Figure 4 is a set of numerical solutions of the time- 
dependent Gross-Pitaeveskii equation that illustrate the ro- 
bustness of the STA for realistic BEC experiments. We con- 
sider a box trap with Gaussian barriers and in all numerical 
simulations interactions are kept constant, i.e. go(t ) — go(0), 
deviating from the ideal prescription in Eq. (16). The top row 
illustrates the dynamics for a quasi-lD BEC. The (one-body) 



fidelity between the resulting state < I ) (t) and the ground state 

P>/ ic <3T l/rhf-rMrl)/ \|2 . 



For 



of the final box <P J gs is J? r=10 = \ (<P(z)\% s )\ z = 0.911. 
smaller values of 7, the fidelity is even higher = 0.999 as 
expected, given that implementation of the exact STA requires 
a smaller tuning of g\. The bottom row shows the dynam- 
ics of a quasi-2D cloud, which requires no interaction tuning 
in a STA, but is more sensitive to the smoothness of the box 
boundaries, J?y = io = 0.988. 

It is noteworthy that in the Thomas-Fermi regime, the ki- 
netic term contribution can be neglected and it is possible to 
induce an exact self-similar dynamics (and a STA) exclusively 
with the help of an external field. Then, the scaling ansatz is a 
solution of the TDGPE provided that 
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This regime is particularly robust against the smooth bound- 
aries of physically realizable box potentials. The simula- 
tions correspond to the most delicate regime with moderate 
mean-field interactions, both far from the non-interacting and 
Thomas-Fermi limits. 
Discussion 

In conclusion, we have presented a method to drive an ultra- 
fast dynamics in a time-dependent box trap which reproduces 
the adiabatic result at the end of the evolution. The method 
is assisted by an auxiliary external harmonic potential which 
provides the speed-up and is applicable to a large family of 
both non-interacting and interacting many-body systems sup- 
porting dynamical scaling laws, where it not only leads to a 
robust expansion of the density but also preserves the non- 
local correlation functions of the initial state, up to an expan- 
sion factor. The proposal is applicable to realistic box po- 
tentials and can be implemented in the laboratory with well- 
established technology. Its applications range over all scenar- 
ios requiring a shortcut to adiabaticity, i.e., probing strongly 
correlated phases, preventing decoherence, the effect of per- 
turbations and atomic losses. The method can be directly ap- 
plied as well to ultrafast population-preserving cooling meth- 
ods, quantum heat engines and refrigerators [ ] providing 
an alternative to the paradigmatic model of a quantum piston 
[14]. 
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